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Abstract 

We study the predictions of the simplest SU(5) grand unified model within 
the framework of minimal supergravity, including constraints from the radia- 
tive breaking of electroweak symmetry. As a consequence of the unification of 
the 5-quark and r-lepton Yukawa couplings, the top quark mass is predicted 
to be close to its fixed point value. We delineate the regions of the super- 
gravity parameter space allowed by constraints from the non-observation of 
proton decay and from the requirement that the LSP does not overdose the 
universe. These constraints lead to a definite pattern of sparticle masses: the 
feature unique to Yukawa unified models is that some of the third generation 
squarks are much lighter than those of the first two generations. Despite the 
fact that all sparticle masses and mixings are determined by just four SUSY 
parameters at the GUT scale (in addition to mt), we find that the signals for 
sparticle production can vary substantially over the allowed parameter space. 
We identify six representative scenarios and study the signals from sparticle 
production at the LHC. We find that by studying the signal in various chan- 
nels, these scenarios may be distinguished from one another, and also from 
usually studied "minimal models" where squarks and sleptons are taken to be 
degenerate. In particular, our studies allow us to infer that some third gen- 
eration squarks are lighter than other squarks — a feature that could provide 
the first direct evidence of supergravity grand unification. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is increasingly being regarded as the most compelling of the 
various extensions of the Standard Model proposed to date. SUSY relates properties of 
bosons and fermions resulting in a new level of synthesis. Supersymmetry also provides a 
technically natural framework for incorporating the Higgs mechanism crucial for the breaking 
of electroweak symmetry. On the experimental front, while the measurements of the gauge 
couplings by the experiments at LEP are inconsistent with minimal models of unification 
of gauge interactions, these measurements are in good agreement with the predictions of 
the simplest supersymmetric grand unified model 0. As an added bonus, supersymmetry 
provides a natural candidate for cosmological dark matter 0. 

Despite these many attractive attributes of supersymmetric theories, no compelling 
model has as yet emerged. This could be due to the fact that the dynamics of supersymme- 
try breaking is not yet understood. In phenomenological analyses, supersymmetry breaking 
is incorporated by the introduction of soft SUSY breaking terms which do not jeopardize 
the stability of the electroweak mass scale in the presence of radiative corrections from very 
massive degrees of freedom not included in the SM (these could be the particles of a grand 
unified theory (GUT), flavour physics at intermediate scales or effects of quantum gravity 
at the Planck scale). Even in the simplest low energy effective theories, the incorporation of 
all possible such soft breaking terms ^ consistent with the known low energy symmetries 
leads to a proliferation of new parameters, greatly limiting the utility of this approach for 
SUSY phenomenology. 

This situation can be greatly ameliorated if the symmetries of the interactions that lead to 
the breaking of supersymmetry can be guessed at. For example, in a SUSY GUT, if SUSY 
breaking is assumed to occur above the scale of GUT breaking, the soft SUSY breaking 
masses and couplings of all particles in the complete multiplet are related. Supergravity 
(SUGRA) models provide a very attractive realization of this idea of utilizing symmetries 
to restrict the supersymmetry breaking interactions. It is assumed that SUSY breaking 
in the effective low energy theory relevant for phenomenology arises due to gravitational 
interactions, which being universal, allow for only a few independent soft-supersymmetry 
breaking parameters (a common scalar mass, a gaugino mass along with SUSY breaking 
trilinear, and possibly, bilinear SUSY breaking scalar coupling) despite our ignorance of 
the details of the dynamics of SUSY breaking. All this, of course, holds at the ultra-high 
scale of SUSY breaking. In order to use the structure of the couplings thus obtained for 
perturbative calculations at the 100 GeV energy scale relevant to experiment today, we must 
sum the large logarithms arising from the disparity of these scales. This is most conveniently 
done using the renormalization group technique |^ whereby the soft breaking masses and 
couplings at the scale of electroweak breaking 0J|] may be obtained. When mixing effects are 
negligible (as is the case for the first two generations of sfermions) or absent (as for gluinos), 
these masses only differ slightly from the physical (pole) masses of the sparticles 0. The 
real beauty of this framework is that these same radiative corrections drive the desired 
breakdown of electroweak symmetry, leaving colour and electromagnetic gauge invariance 
intact. As a bonus, supergravity models lead to (approximate) degeneracies of the first two 
generations of squarks which is required |T0[ by the non-observation of flavour changing 
neutral currents in the K meson sector. It is interesting that the masses of t- and b- squarks 
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may be significantly different; this can lead to interesting phenomenological consequences 
as we will see below. 

While most phenomenological analyses of supersymmetry embody some of the ideas of 
supergravity GUTs, the observation that the LEP data are consistent with the predictions 
of SUSY SU(5) has led several groups to systematically study the implications of the 
simplest SUGRA SU(5) model with radiative symmetry breaking constraints incorporated. 
It has been pointed out |]TT| that these same data vitiate one of the earliest successes of 
non-SUSY SU(5) GUT — the prediction of the ratio rrib/m^ ~ 3 from b-r Yukawa coupling 
unification. In contrast, the data are in agreement with the Yukawa coupling unification in 
the SUGRA SU(5) framework provided that the top quark mass is close to its infrared fixed 



point [|12H l 



The purpose of this paper is to study the implications of the simplest SU(5) supergrav- 
ity model with radiative symmetry breaking constraints incorporated. We will mainly focus 
on the detectability of the sparticle signals at the Large Hadron Collider (LHC) and espe- 
cially, on whether these signals serve to distinguish this framework from the usually studied 



"minimal SUSY" framework [^,|T^. Here, by simplest, we mean that we assume the fewest 
number of particles and the smallest number of soft-SUSY breaking parameters at the scale 
Mx of GUT symmetry breaking. More specifically, we assume: 

1. A minimal set of soft SUSY breaking parameters at the GUT scale Mx — 2-10^^ GeV, 
i.e. a common scalar mass mo, a common gaugino mass mi/2 and common trilinear 
{Aq) and bilinear {Bq) soft breaking parameters. We do not assume any relation 
between Aq and Bq, though. 

2. Minimal particle content at and below the scale Mx- 

We recognize that these assumptions are sensitive to the unknown physics at high energy 
scales, and may well prove to be incorrect; for instance, the assumption of a common gaugino 



mass might be invalid if there are non- minimal kinetic terms for the gauge fields Our 
point, however, is that at present even this very restrictive model can be claimed to be 
realistic in the sense it is compatible with all experimental constraints both from particle 
physics as well as from cosmology. Moreover, as we will see, the fact that the top quark 
Yukawa couplings are driven to their fixed point values results in supercollider signals that 



are characteristically different from previous studies Despite the fact that our analysis 
entails an extrapolation of physics over a very large range of energy, we believe that the 
simplicity of the model and the characteristic predictions of the scenario provide sufficient 
motivation for our analysis. 

In the next section, we discuss how the two assumptions together with various observa- 
tional constraints from collider and non-accelerator experiments as well as from cosmological 
considerations further constrain the model parameters in this already restrictive framework. 
We also recapitulate how the sparticle spectrum and the SUSY parameters of the low energy 
theory are obtained and present representative examples of these. In section III, we study 
the general features of the various scenarios obtained with our assumptions, and discuss how 
these differ from those usually considered in phenomenological analyses within the Minimal 
Supersymmetric Standard Model (MSSM) framework. In the next section we briefly de- 



scribe our simulation of sparticle production at hadron colliders using IS A JET 7.07 |T8 



and discuss various improvements that we have recently incorporated into the program. In 
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Section V, we discuss various SUSY signals that might be expected at the Large Hadron Col- 
hder (LHC) in the representative scenarios introduced above. We show that an observable 
signal may be expected in various event topologies including in single lepton and dilepton 
+ If J, samples, which are usually thought to be obscured by large SM backgrounds. Here, 
and in the next section, we discuss how the simultaneous observation of a signal in several 
channels can serve to differentiate the SUGRA SU(5) from the MSSM framework, and also, 
to distinguish between the various illustrative supergravity scenarios. We conclude with a 
summary of our results and some general remarks. 

II. THE SUPERGRAVITY SU(5) FRAMEWORK: DESCRIPTION OF 

CONSTRAINTS 

In this section, we first overview the relations between mass parameters at the weak 
scale in the minimal SUGRA model. We then describe various theoretical and experimental 
constraints that we impose, and show that despite these, it is possible to find viable solutions 
even in our theoretically restricted scenario. 

At the high scale Mx-, we start ||^,|^ with the common gaugino mass mi/2, the common 
scalar mass mo, the higgsino-mixing parameter /xq, the soft SUSY breaking trilinear (bilinear) 
couplings Aq (Bq), together with the Yukawa couplings ht, hb and (with hb = h^. at scale 
Mx] see below). The absence of any intermediate scales implies that all parameters at the 
weak scale can be determined uniquely from the input parameters at Mx, by solving a set of 
renormalization group equations (RGE) 0. Of course, the above set of parameters should 
reproduce the correct size of symmetry breaking, or Mz, rrit (as a input), rrib and at the 
weak scale; this leads to constraints between the parameters at Mx- We do not assume any 
relationship between Aq and Bq. In practice we use the ratio of vacuum expectation values 
tan (3 as an input parameter. Using the equations describing the minimization of the Higgs 
potential at the weak scale then allows us to express B and at the weak scale in terms of 
mz and tan/? 0; note that the renormalization group equations for B and /i decouple from 
all other RGE. The constraint of proper electroweak symmetry breaking therefore reduces 
the number of free input parameters to five, which we take to be mg, mi/2, ht and 
tan/3; in addition, the sign of needs to be specified. 

These considerations lead to simple relations between gaugino mass parameters (Mj) and 
first and second generation sfermion masses at the weak scale. The RG running of the masses 
of strongly interacting particles implies large positive corrections to squark and gluino mass 
parameters, while corrections to the electroweak interacting particles stay relatively small. 
Numerically, m? ~ rriQ + GmJ^g; the soft SUSY-breaking gaugino masses are related by 
Ml ~ 0.4mi/2, M2 ~ 0.8mi/2 and M3 ~ 2.7mi/2. 

The Higgs mass parameters and third generation squark masses receive significant ra- 
diative corrections from Yukawa interactions in addition to those from gauge interactions. 
At Mx our ansatz requires the two Higgs mass parameters rriHi and 171^2 to be equal: 
m'jj_ = ml + /Xq. As is well known, this equality implies that the tree-level potential can- 
not simultaneously lead to nonzero vacuum expectation values and be bounded from below. 
Quantum effects, i.e. the running of the mass parameters of the scalar potential, thus be- 
come crucial for gauge symmetry breaking. Radiative SU{2) x f/(l) symmetry breaking is 
in fact a natural prediction of this model, given the large size of the top Yukawa coupling. 
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More explicitly, receives a negative radiative correction proportional to hf times some 
function of soft SUSY breaking parameters when going towards the weak scale, due to its 
interaction with the top quark. On the other hand m|^^ is insensitive to ht, thereby lifting 
the equality between Higgs mass parameters. If hf is very large, m^^^ becomes negative well 
above the weak scale, resulting in too large a value of mz or even an unbounded Higgs 
potential. This can be compensated by choosing to be large. In other words, has to 
be increased along with ht in order to maintain Mz at its correct value. 

We now turn to the constraints coming from our assumption of a minimal SUSY SU{5) 
GUT. The minimality of the SU{5) Higgs sector immediately implies the equality of the 
bottom and tau Yukawa couplings at the unification scale. Several recent studies [|T^-[T^ 
have shown that this is only compatible with the measured values of b quark and r lepton 
masses if the top Yukawa coupling is large, i.e. close to its "fixed point" value |]19[. This 



conclusion persists even if the ratio ht/hr is varied about unity by 10% to account for 
uncertainties due to threshold effects. This immediately leads to the prediction 

mtimt) ~ 190 GeV ■ sin /?, (2.1) 

where mt{mt) is the running (MS) mass. Since we only use 1-loop ROE everywhere, we 
simply interpreted eq. ( ^.1| ) to mean 

htiMx) = 2. (2.2) 



Taking the on-shell top quark mass rrit < 190 GeV [Q, we find tan/3 < 2.5 |T2[, in general 
agreement with the constraints on it from nucleon decay [^] which we now turn to. 

In SUSY GUTs, the exchange of higgsino triplets can lead to a dangerously short lifetime 
of the proton. In minimal SUSY SU{5) its exact value depends on many parameters. In 
particular, the lifetime increases with increasing mass of the higgsino triplet. Since the 
dimension-5 operators induced by higgsino triplet exchange need to be "dressed" by loops 
involving squarks and gauginos (chargino loops usually being dominant), contributions from 
(s) quarks of different generations in the loop add coherently. This opens the possibility 
to increase the total lifetime, at the cost of some fine-tuning, by cancelling contributions 
involving third generation (s) quarks against those involving only the first two generations. 
The "dressing loop function" can be decreased, and hence proton decay can be suppressed, by 
choosing scalars (squarks and sleptons) to be considerably heavier than gauginos. For given 
values of the other parameters the lower bound on the proton lifetime therefore translates 
into an upper bound on the ratio rni/2/ml. In view of the various uncertainties p2[ of this 
calculation we have conservatively required ||21|| , 



mo > min(300 GeV, Smi/s). (2.3) 

The coefficient 3 on the right hand side has been chosen so that the contribution rriQ to m? 
described earlier dominates the one oc which implies rriq S> mi/2. 

Finally, we require our model to be cosmologically acceptable. Because the MSSM has a 
conserved R parity, the lightest supersymmetric particle (LSP) is stable. Thus the LSP can 
be a good candidate for the missing cosmological mass density, or Dark Matter (DM). Large 

means that the lightest neutralino, which is also LSP in this scenario, is dominantly a 
gaugino. The constraint eq. (|2.3|) implies that LSP annihilation via the t— channel exchange 



5 



of sfermions is suppressed, leading to a dangerously large relic abundance (or small age 
of the Universe) unless LSP annihilation is enhanced by the proximity of an s— channel 
pole. This is not as unnatural as it sounds, since recent studies have shown very large 
pole enhancements even if the LSP mass is several GeV below half the mass of the boson 
exchanged in the s— channel. Our large values of fi and mo and the rather small value 
of tan/3 force most Higgs bosons to be heavy, the exception, of course, being the lighter 
neutral scalar Hi. This is because m|^^ does not receive large negative radiative correction 
from Yukawa interactions if tan (3 is small, while = fn\j^+m'\j^ with Tn\j^ > —m\/2 from 
the minimization of the Higgs potential. The condition eg. ( f^.3[ ) then implies > lOmLSP- 
Therefore only the Z boson and the light Higgs boson can contribute significantly to LSP 
annihilation. 

In our computation of the Higgs sector, we have included one loop corrections from top 
and bottom (the latter are small) Yukawa interactions to both masses and trilinear Higgs 
interactions using the effective potential 1^^; note that these corrections are not contained 
in the RGE discussed above. Despite the large top Yukawa coupling, we find that tjih^ is 
bounded from above by about 110 GeV for the values of tan jS we are considering. Moreover, 
Hi is always SM-like in this scenario, due to the large masses of the other Higgs bosons. 
Hence the LEP bound of 62.5 GeV applies for it as well. This then leads to bounds on the 
mass of a cosmologically acceptable LSP, and hence on the gaugino mass: 

60 GeV < |mi/2| < 130 GeV. (2.4) 

Finally, we impose the limit 

mo < 1 TeV (2.5) 
in order to avoid excessive fine tuning. 



We should mention a technical complication that arises from the constraint (2^) as a 
result of which squarks are substantially heavier than gluinos. It has been shown that the 
pole gluino mass receives additional corrections proportional to \og{mq/mg). Since one of 
the features of these scenarios, as we will see, is a rather light gluino, these extra corrections 
- which can be as large as 30% - can be significant. We have incorporated these into our 
program using the formulae in Ref. 0. 

Our apparently innocuous assumption of minimal particle content at the GUT scale led 
to the various constraints on model parameters as we have just discussed. Most analyses 
require either the equality of b and r Yukawa couplings at Mx, or an acceptably long lifetime 
for the proton, but not both. It seems to us that these two constraints are really on very 
similar theoretical footing in that they can both be weakened or even avoided completely if 
one allowed a more complicated Higgs sector and/or a larger GUT group ||2^. On the other 
hand, we do not try to find correlations between the value of the strong coupling constant 
and the overall size of SUSY breaking masses. Given the present (mostly theoretical) error 
of as{Mz), and in view of further uncertainties at the high scale we do not consider any 
"bounds" on sparticle masses derived in this fashion to be reliable at present. We now turn 
to a discussion of the characteristics of the fixed point scenario outlined above, and give 
some illustrations of the spectra and the sparticle mixing patterns that may be anticipated 
in this class of models. 
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III. THE SUPERGRAVITY SU(5) FRAMEWORK: SPARTICLE MASSES AND 

MIXING PATTERNS 



The characteristic feature of the fixed point scenario of the last section is the large top 



Yukawa coupling, eq. (|2.1| ). This, as we have seen, leads to > niQ, assuming the radiative 
breaking mechanism of electroweak symmetry. Together with eq. ( p.3|) , this then implies 
that the lighter (heavier) charginos and neutralinos are gaugino- (higgsino-) like in character. 
The limit eq. ( p.4| ) implies that gluinos cannot be much heavier than about 400 GeV and 
may even be light enough to be discovered at the Tevatron. We have also seen that the 
pseudoscalar (Hp), the charged (-ff^) and the heavier neutral scalar (Hh) Higgs bosons all 



tend to be heavy, so that these will not be produced via cascade decays of gluinos p6 
Finally, we make the obvious remark, rriq > rrig so that the squarks of the first two families 
dominantly decay to gluinos. 

Another important consequence of the large top family Yukawa interactions is the sizable 
reduction of the SUSY-breaking masses of scalar top quarks, and further, the mass squared 
of is reduced twice as much as that of t^. By SU{2) invariance this means that the 
mass of left-handed scalar bottom quarks will also be reduced compared to squarks of the 
first two generations. The top squark masses are further altered by squark mixing, which 
depresses the lighter state even further. This non-degeneracy of third generation squarks, 
which is one of the characteristic features of minimal SUGRA models with radiative gauge 
symmetry breaking is accentuated by the large value of ht in the model. 

In spite of the rather strong constraints on parameter space, we found that there is 
still considerable room for variation in the expected phenomenology. While condition ( p.3|) 
implies that masses of the first and second generation squarks are well above the gluino mass 
no such statement can be made for the lighter stop eigenstate ti. In particular, it is possible 
that the two body decay g —>■ t + ti is allowed |^, in which case it will have a branching 
ratio near 100%, being the only two-body decay mode of g. We remind the reader that 



it has been shown in Ref. |^ that signals from a relatively light t-squark may well have 
escaped detection at the Fermilab Tevatron, although with the accumulation of a large data 
sample, its detection may be possible . 

Despite the various constraints that we have imposed it is impossible to do an exhaustive 
scan of parameter space. Moreover, for the purposes of this study, the main question is not 
whether it is possible to detect supersymmetry at the LHC (for the light spectrum we have 
this has already been answered affirmatively |16|), but whether it is possible to distinguish 



the present class of models from the usually studied MSSM, and further, whether it is possible 
to distinguish the various scenarios (discussed shortly) from one another. Toward this end, 
we have attempted to choose model parameters (respecting the constraints discussed in Sec. 
II) that accentuate the qualitative differences in the signals. 

At hadron colliders, squarks and gluinos are the most copiously produced particles. Since, 
within our framework, the first two generations of squarks dominantly decay to gluinos, we 
expect that the gluino decay patterns fix the broad features of the phenomenology. This 
immediately leads us to consider two classes of scenarios: (A) where the decay g ^ t + ti 
is allowed, and (B) where it is forbidden. Case A naturally subdivides into two cases 
determined by the decay pattern of the lighter t-squark. In case Al, we assume that the 
decay ti — > bWi is accessible; since this is the only tree-level two body decay, its branching 
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fraction is essentially 100%. In case A2, we assume that the decay ti bWi is kinematically 
inaccessible; in this case the stop decays [^] via one loop diagrams into a charm jet and an 
LSP: ti ^ cZi. 

Turning to the heavier stop, case B, we have attempted to diversify the phenomenology 
by choosing combinations of mo and mi/2 at their maximum and minimum values discussed 
in Sec. II. This leads to four cases B1-B4. 

• Bl: mi/2 and mo small, 

• B2: mi/2 small; mo large, 

• B3: mi/2 large; mo small, 

• B4: mi/2 and mo large. 

Notice that for the case Bl where both squarks and gluinos are at their minimum phe- 
nomenologically acceptable values, the gluino may even be accessible at the Tevatron. Our 
point in analysing this case is more to illustrate the diversity of phenomena that could re- 
sult even in this restricted framework, and to study the range of parameters over which 
distinguishable signals may be possible at the LHC 

The input parameters, the resulting spectra and the values of fi at the weak scale in 
these six cases are shown in Table I. Several features are worth pointing out. 

• We see that the gluinos are relatively light. As noted above, gluino signals in case Bl 
(and, possibly, B2) may even be observable at the Tevatron. 

• It is striking to see that in case A2, m^^ is only 84 GeV; a significant number of ti-pair 
events should, therefore, be present in the existing data sample at the Tevatron. As 
pointed out in Ref. such a light stop squark may well have escaped detection at 



the Tevatron especially since the LSP is rather heavy in this case. It should, however, 
be remembered that both the experiments at the Tevatron have since accumulated 
at least 15-20 pb~^, so that the detectability of the stop at the Tevatron merits a 
reevaluation |2^. It would indeed be interesting to study whether stop squark signals 



in cases A and B2 might be detectable at the Tevatron. 

The lighter chargino and the second lightest neutralino in Bl and B2 should be kine- 
matically accessible at LEP II. For the other cases in Table I, their accessibility will 
obviously depend on the energy that LEP II is finally able to attain. Since, as remarked 
above, both Zi and Z2 are dominantly gaugino-like (and selectrons are heavy), their 
production at LEP II will be suppressed. For cases Bl and B2, the trilepton signal 
from W1Z2 production at the Tevatron should be present at observable rates. In the 
other cases illustrated, the trilepton cross section (after experimental cuts) is ~1 fb 
pl[] , and hence, is unlikely to be observable. This is because we have m^ ~ m^ so 



that, unlike as in no-scale models, the leptonic branching ratios are not dynamically 
enhanced. 

The light Higgs boson is essentially the SM Higgs boson. It should, therefore, 
certainly be detectable at LEP II in case Bl, while the detectability in the other cases 
will be sensitive to the machine energy that is ultimately achieved. 
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We note that the top quark is necessarily within the reach of the Tevatron, and further, 
its decays are as given by the SM except in case A2 where the decay t ^ tiZi has a 
branching fraction of < 10%. 

Finally, we see that as noted above, m^^ is indeed considerably smaller than the mass 
of the first two generations of squarks. This implies that three body gluino decays 
mediated by bi exchange are enhanced relative to those mediated by other squarks. 
As a result, gluino events in cases B may be expected to have a higher multiplicity 
of B mesons than would be expected if all squarks are taken to be degenerate as in 
most MSSM studies. It is amusing to note that gg events may be even richer in 5's 
for cases A, since then, decays of t quarks from g — > tit as well as the decays ti —>■ bWi 
(in case Al) lead to B meson production. We also note that because 6l is relatively 
light, Z2 bbZi decays can also be considerably enhanced relative to Z2 ddZi 
decays. We thus see that the presence of light third generation squarks leads to SUSY 
events that should be rich in B mesons. Such events may also occur within the MSSM 
framework where cascade decays of gluinos to Higgs bosons can result in events with a 
high multiplicity of b quarks [^. However, these cascade decays can occur only when 
m^^ — m^^ exceeds the lighter Higgs boson mass which, for the SM Higgs bosons, is 



bounded below by 63.5 GeV ||33|. Thus, within the MSSM framework, these events 
occur only when gluinos are rather heavy. In other words, the detection of 5-rich 
gluino events with event kinematics corresponding to gluino masses smaller than 400 
GeV could signal a light third generation of squarks. 



Next, we briefly describe the improvements that we have made in ISAJET 7.07 |18 
to simulate supersymmetry events at hadron colliders that we use to compute event rates, 
topologies and kinematic distributions for the six cases listed in Table I. 



IV. SIMULATION OF SUPERSYMMETRY AT HADRON COLLIDERS 

The main production processes and decay modes of all sparticles as given by the minimal 
supersymmetric model have been incorporated into ISAJET [jl8[ which we use for all our 
simulations. In order to allow for the possibility that the third generation of squarks may be 
substantially lighter than the first two generations, which is the distinguishing characteristic 
of the scenario, we have had to make several modifications to the code. In ISAJET 7.07, 

1. we have included the production and decays of t-squarks. 

2. We have modified the input parameters to the program. In addition to the SM param- 
eters (including nit) the current low energy inputs in terms of which sparticle masses 
and decay patterns are computed are (m^, nig, rrii^, rn^^, mp, tan/?, /x, niHp, ^i^, 
TTii^, At, ^1^1 Ah). The last five entries are new modifications and allow for indepen- 
dent mass and mixing pattern for the third squark generation. The stop and sbottom 
masses are the respective soft SUSY-breaking mass parameters and not the physical 
squark masses. Since b^ and belong to the same SU(2) doublet, m^^ = m^-^. We 
have neglected D-term contributions that break the degeneracy between the L- and R- 
type squarks of the first two generations but have retained these in the computation of 
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slepton masses — this is the only reason for the difference in the and vij^ masses. For 
heavy squarks, this D-term sphtting is generally small, and so, makes little difference 
in the computations of squark production cross sections or gluino, chargino and neu- 
tralino decays mediated by virtual squarks. This splitting can be very important when 
squarks and gluinos are close in mass, since then some of the two-body gluino decay 
modes may become accessible. In the present study, this is never the case because of 
the constraint (12. 31). 



3. The input gluino mass is the physical pole mass as opposed to the running MS 
mass that enters the gaugino-Higgsino mass matrices. In all the cases in table I, the 
pole mass is somewhat larger than the MS mass, primarily because the squarks are 
heavier than gluinos. For Case B2, where the gluino is light and squarks are at their 
maximal value, this difference is almost 50 GeV. 

4. In the computation of the Higgs boson masses and couplings, we have, as already 
mentioned, incorporated radiative corrections from both top and bottom Yukawa in- 
teractions, including mixing between third generation squarks. We have not included 



effects of virtual gauge boson and gaugino loops, but these are known to be small [34 



5. We have supplemented the various sparticle decay formulae to account for non- 
degenerate third generation sparticles. This can be especially important for gluino 
three body decays, as well as the corresponding neutralino decays. We have also in- 
cluded contributions to three body Z2 decays mediated by Higgs bosons; these can be 
important for large values of /i and niq for which the contributions from both Z and 
squark mediated decays are dynamically suppressed. As a result, Z2 decays to heavy 
quarks and leptons may be enhanced. 

For every set of input parameters, the routine ISASUSY calculates all sparticle masses 
along with branching fractions for all their allowed decay modes. Of course, the input 
parameters are correlated by the minimal SU(5) supergravity constraints discussed above. 
We then use ISAJET to generate supersymmetric particle pairs in pp collisions. In each 
run, ISAJET is set up to simultaneously generate all allowed sparticle pairs in numbers 
proportional to the corresponding pair production cross section. We use the EHLQ Set I 
structure functions |^ for our computations of cross sections at the LHC. The produced 
sparticles are then decayed via various cascades with branching fractions as obtained from 
ISASUSY into final states involving quarks, leptons, gluons and photons. Radiation off 
initial and final state partons as well as heavy flavour fragmentation is included in ISAJET. 
Daughter quarks and gluons are hadronized, and unstable particles decayed into the final 
state. Finally, underlying soft scattering and hadronization of beam remnants is modelled 
in our simulation. 

For each of the six cases in Table I, we use ISAJET to generate a sample of 50K SUSY 
events by pp collisions at i/i = 14 TeV. The relative contributions of the different SUSY 
processes to this event sample is shown in Table II. We see that for the light stop cases 
Al and A2, the production of stops dominates, and furthermore, that the stops have a 
significant production rate even in most of the cases B. The production of other flavours 
of squarks is suppressed relative to that of gluinos primarily because of Eq.( p.3| ). We have 
also shown the branching fractions for the main decay modes of sparticles that are most 
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abundantly produced (either directly or via cascades) at the LHC In cases Al and A2, 
the gluino essentially always decays to t-squarks whereas in cases B, Wi and Z2 decays 
are favoured We also see that gluino decays to b quarks are enhanced relative to the 



corresponding decays to d quarks. This is because 6^ is lighter than the other down type 
squarks. This enhancement is small for cases Bl and B2 because the Zibhi coupling is 
suppressed by accidental cancellations — these cancellations do not occur for gluino decays 
to Z2. Notice also that the decays Z2 —>■ bbZi are also enhanced relative to Z2 — »• ddZi 
decays, in part for the same reason, but also because of exchange diagrams. Finally, we 
have also shown the branching fraction for the leptonic decays of Wi and Z2 which enter 
the computation of the multilepton signals. Notice that the neutralino branching fraction 
exhibits substantial variation while that for the chargino is essentially determined by the W 
branching fraction except in Case A2 where the chargino decays via the two body tib mode. 

In order to simulate experimental conditions at the LHC, we impose cuts described in 
the next section and then classify the signals according to their lepton content as in Ref. 
T6| . Since the efficiency for SUSY events to pass these cuts depends strongly on the masses 



and decay patterns of the produced sparticles, the contribution of the various processes after 
cuts is, in general, quite different from the relative rates in Table II. To compare the features 
of sparticle production in the supergravity SU(5) model under study, with the corresponding 



production in "minimal SUSY" models usually p!6| , p!5| used to study SUSY signatures, we 
repeat our simulation for four illustrative cases. Here, we take all squarks and sleptons to 
be degenerate, and fix the gluino mass at 300 GeV, a value similar to rrig in the six cases in 
Table I. Since the signals are known to be sensitive to mq/rrig, we illustrate our results for 
rriq = TUg + 20 GeV and rriq = 2mg. To vary the electroweak gaugino content we study two 
values of fi. Specifically, the parameters for these four cases are: 

• BTWl: m~g = 300 GeV, = 320 GeV, = -150 GeV, 

• BTW2: m-g = 300 GeV, rrig = 600 GeV, /i = -150 GeV, 

• BTW3: m-g = 300 GeV, rrig = 320 GeV, fx = -500 GeV, 

• BTW4: m-g = 300 GeV, m^- = 600 GeV, fx = -500 GeV. 



For comparison purposes, and in keeping with parameters of Ref. [0, we have fixed 
tan (3 = 2, again a value similar to the cases in Table I, and chosen mt = 140 GeV, mHp = 500 
GeV, with the A-parameters set to zero. A discussion of our simulation along with the 
various results forms the subject of the next section. 



V. SUPERGRAVITY SU(5) SIGNALS AT THE LHC 



As emphasized in Ref. several interesting signals should simultaneously be present 
if any new physics signal at the LHC is to be attributed to squark or gluino production. 
These signals include the canonical I^j, events, the like-sign dilepton events and multilepton 
events. These leptonic events should contain substantial jet activity. The production of 
charginos or neutralinos in association with a squark or gluino yields the same event topolo- 
gies, while chargino-neutralino production could, in addition, lead to multilepton events free 
of jet activity. In this section, we study the rates and distributions for these various event 
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topologies for the six supergravity cases in Table I and compare these with the "minimal 
super symmetry" test cases listed above. 

We use the toy calorimeter simulation package ISAPLT to model detector effects. We 
simulate calorimetry with cell size Arj x A0 = 0.1 x 0.1, which extends between —5 < ?7 < 5 
in pseudorapidity. We take electromagnetic energy resolution to be 10%/a/Et © 0.01, 
while hadronic resolution is 50%/^^ © 0.03 for \r]\ < 3, and 100%/^^ © 0.07 for 
3 < \ri\ < 5, where © denotes addition in quadrature. Jets are coalesced within cones 
of i? = \/ Arj'^ + A(j)'^ = 0.7 using the ISAJET routine GETJET. Hadronic clusters with 
Ej- > 50 GeV are labelled as jets. Muons and electrons are classified as isolated if they 
have pt > 20 GeV, \ri{i)\ < 2.5, and the visible activity within a cone of R = 0.3 about the 
lepton direction is less than i?T(cone) = 5 GeV. Finally, we assume (perhaps optimistically) 
that it is possible to identify 6-jets with px >20 GeV within |?7| < 2 with an efficiency of 
40% via micro-vertex tagging. 

Following Ref. |1^, we group the various signal categories according to the number of 
isolated leptons, and further, require I^j. > 100 GeV for the SUSY signal. The various 
signals are, 

(z) the I^rp signal consisting of Oi + rij jets + I^j, events. We require I^j. > 150 GeV, 
Uj > 4 and the transverse sphericity St > 0.2; 

(a) the single lepton signal consisting of events with one and only one isolated lepton 

(1£ + + 

(m) events with two isolated leptons {2i + jets + I^j,), which can subdivide into opposite 
sign dileptons (OS) and same sign dileptons (SS). 

(iv) events with three isolated leptons {3i + jets + ^y); 

(v) and events with four isolated leptons (4£ + jets + J^rp). 

Table III gives an overview of the cross sections for these event topologies for the six 
SUGRA cases in Table I along with the cross sections for the major SM backgrounds. Also 
shown are the corresponding cross sections for the four "MSSM" cases introduced in Sec. 
IV. The cross sections shown here include the various cuts mentioned above, as well as 
additional cuts (discussed below) for the the li and the OS dilepton signals which are 
known to have large SM backgrounds. As we will see, these additional cuts are essential to 
suppress the SM background, particularly to the li and £~^£~ signals to acceptable levels. 
We see from Table III that for all but the li signal in Case B4, the signal to background 
ratio exceeds unity, and is frequently considerably larger. Furthermore, we see that the 
cross sections in the 1£ and dilepton channels range from 1-100 pb so that 10^-10^ signal 
events may be expected annually at the LHC even for a luminosity of 10^ pb~^/year. We 
also see that except for case A2, several thousand trilepton events may be expected every 
year at the LHC. We have been unable to identify any significant SM physics background 
to this signal. Finally, we note that while there is some variation of the signals for the 
different scenarios in Table HI, these cross sections alone do not enable us to unambiguously 
distinguish between all the various cases, or for that matter even between the SUGRA and 
the "MSSM" cases. This requires a more detailed study of the various signals which we now 
turn to. 
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A. Ifij. signal 



This is generally regarded as the canonical signal for supersymmetry. Backgrounds in 
this channel can be very dependent on detector details, especially forward calorimetry, and 
realistic BG estimates require a detailed calorimeter simulation. Nevertheless, we estimate 
backgrounds from several major SM sources using our toy calorimeter and ISAJET: tt pair 
production, for rrit = 160 GeV, W +jets production, Z + jets production, and heavy flavor 
{hh and cc) production. Considerable background rejection with relatively little cost to signal 
can be gained by requiring, in addition, 

• SIE't-I > 700 GeV, where the sum is over all calorimeter cells. 

The impact of this cut on the signal and the various backgrounds is shown in the Ifj, 
column of Table IV. The values (inside) outside the parenthesis indicate the signal cross 
section (before) after this additional cut. While this cut is not essential in that the signal to 
background ratio substantially exceeds unity even before the cut, we see that implementing 
the cut results in a significant improvement of this ratio with only a modest loss of signal. 
The largest I^j. cross sections come from the low mi/2 (light gluino) cases, where the total 
production cross sections are largest. In order to facilitate distinction between the different 
scenarios in Table III, we exhibit various characteristics of the events in Table V. In 
columns 3-5, we list the fraction of events with jet multiplicities rij = 4 — 5, rij = 6 — 7 
and rij > 8. The most noticeable difference occurs in the nj > 8 channel, where we see 
several SUGRA cases yield large fractions of events with very high jet multiplicity. This is 
especially true of cases Al, A2 and B4. For instance, in Al, many of the events come from 
gluino pair production, where g ^ tii ^ bqq' + bqq'Zi, which can yield up to twelve 'parton 
level' jets. Of course, it is very difficult for the SM background processes to yield such high 
jet multiplicities. 

We also list in Table V the fraction of events containing identifiable displaced 6-vertices, 
Ufy. An important point is that supersymmetry events often have high multiplicity of B 
hadrons ITOJH^ in the central region. This is even more so in the SUGRA cases we have 



generated. We see that in the SUGRA cases, from 25% — 65% of the events have at least 
one identifiable B vertex, compared to 15% — 25% for the BTW cases (of course, larger B 
multiplicities are also possible in the BTW cases when cascade decays of gluinos into Higgs 
bosons become operative. However, this occurs for considerably larger values of rrig.). Also, 
the SUGRA cases can have large fractions of events with Uh > 2 and Uh > 3. There are 
several reasons for this: 

• in case Al with the richest B content, each gluino in gluino pair events decays via 
g tti, which automatically yields at least four i?'s per event since the t-quark as 
well as the ii squark decay into a b, 

• the fact that 6^ is generally lighter than the first two generation squarks means that 
three body gluino decay via a virtual 6^ is enhanced, yielding a large branching fraction 
for -g^bbZ, i, 

• neutralino three body decays such as Z2 bbZi are frequently enhanced, either be- 
cause of the smaller bi mass, or due to the importance of Higgs boson mediated 
diagrams discussed in Sec. IV. 
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Of course, the exact b niTiltiplicity will depend on the tagging efficiency that is ultimately 
attained. The ability to tag B vertices may also be very helpful in obtaining an enriched 
SUSY sample, and thus aid in SUSY particle mass estimates. 

We also show in Table V the mean summed scalar Et and mean in I^j, + jets events. 
Use has already been made of T,\Et\ to enhance SUSY signal against background. The 
scalar Et distribution is somewhat sensitive to the magnitude of the SUSY particle masses 
produced in the subprocesses although this is not clear from just the mean values listed 
in Table V. We also remark that the differences between the various cases are somewhat 
reduced by the hard T,\Et\ cut. We show in Fig. 1 the px distribution of the fastest jet in 
^j, events, for four SUGRA cases, and the sum of backgrounds. The SUGRA cases have px 
which peaks around ~ 200 — 250 GeV, with a broad tail extending to higher pt's, unlike SM 
background. Since squarks of the same mass are expected to lead to harder jets than gluinos, 
it is natural to ask whether squark production shows up as a shoulder in this distribution. 
A glance at Fig. 1 shows that this will be difficult. We remark though that the hard cuts 
that we have imposed to enhance the signal selects out events with similar kinematics and 
so reduces the differences between the various signals. 

B. Single isolated lepton events 

Little effort has gone into searches for SUSY in the single lepton channel, due in part 
to large backgrounds expected from W + jets production and ti production. Within our 
framework, the production rate of (the rather light) gluinos is truly enormous at the LHC 
so that one may hope to be able to make rather stringent cuts to eliminate the backgrounds, 
leaving behind an observable signal. A glance at the cross sections in parenthesis listed in 
the li column of Table IV shows that the single lepton cross section is smaller than the SM 
background except in case Bl. To facilitate further separation between the SUSY signal 
and the SM background, we show in Fig. 2 a) the and b) transverse mass mT{i, Ifx) 

distributions for four SUGRA cases Al, A2, Bl and B4, along with the summed background 
distribution. The scalar Ej- distribution is substantially harder than SM background. Since 
both the top and the W backgrounds dominantly contain one leptonically decaying W boson, 
the transverse mass distribution for the background rises to a Jacobian peak at rriT = Mw 
while that for the signal is relatively smooth. Hence, to maximize signal against background, 
we further require: 

• S|^t| > 700 GcV, and 

• mT{£, $t) < 60 GeV, or mrii, $t) > 100 GeV. 

The resultant cross sections are shown in Tables III and IV, in the li column. The total SM 
background is ~ 12 pb, while signal ranges between 10 — 77 pb. In order to get some idea of 
the hardness of the leptons in the li sample, the Pt(^) distribution after all the above cuts 
is shown in Fig. 2c. In case Bl, where single leptons are apt to come from light chargino 
decay, a distortion of the low prii) distribution can occur; for other cases, where leptons 
come from t-quark decay, the distribution more nearly tracks the shape of the background. 

Various characteristics of single lepton events are detailed in Table VI. The general 
features are similar to those of Table V. The jet multiplicity is again large, especially for 
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events with larger gluino masses, while jet multiplicities for SM backgrounds arc substantially 
smaller. Furthermore, the SUSY signal events, especially for cases Al, A2 and B4. are 
characterized by large fractions with visible displaced B vertices. The SUGRA cases, except 
for case B2, all have > 2 substantially larger than the corresponding BTW values, and 
hence could yield evidence for a SUGRA-type of sparticle spectrum, where some of the third 
generation squarks are far lighter than squarks of the first two generations. 

C. Opposite-sign (OS) isolated dilepton events 

As for the case of single lepton events, searches for SUSY in the isolated opposite sign 
dilepton channel have rarely been considered due to a large expected background from tt 
events. Within our framework, however, we see from the cross sections in parenthesis in 
the OS column in Table IV that, with just the canonical cuts discussed above, the dilepton 
signal exceeds that background in the favourable cases Al, Bl and B2. Furthermore, it is 
only between a factor 2-6 smaller than the background in the other cases. In Fig. 3a we 
plot the SlE'rl distribution, for which the signal gives a substantially harder spectrum than 
background. We see that the additional requirement, 

• El^rl > 700 GeV, 

reduces the background by a factor of seven, with only modest loss of signal, particularly for 
the cases A2, B3 and B4 where the signal was considerably below background. The cross 
sections after all cuts are shown in Tables III and IV, in the OS columns. Signal exceeds 
background for cases Al, Bl and B2, and is comparable in the other SUGRA cases. The 
Pt{^) distribution for fast and slow Icptons after the above cuts is shown in Fig. 3b and 
3 c. The overall hardness of these distributions depends in part on the SUSY sources of 
the dileptons. In cases Al and A2, many of the leptons are quite hard, coming from top 
quark decay, while in cases B2 and B4 the leptons frequently come from three body chargino 
and neutralino decay, and exhibit a softer spectrum. For OS dileptons produced in gluino 
and squark events, cither the leptons each separately come from the two SUSY particles 
produced in the hard scattering, or the two leptons can come from a single SUSY particle 
cascade, such as ^ — > — > iiZi. In the former case, the dilepton pair is expected to have 
a large opening angle in the transverse plane, whereas in the latter case, the opening angle 
is expected to be small. This is borne out in the behavior of Fig. 4. In cases B2 and 
to a smaller extent in case B4 (cases Bl and B3 are similar to Bl and B4, respectively), 
where there is substantial neutralino production via cascade decays, the dileptons come 
with a small opening angle. In case Al the dileptons are somewhat peaked at large opening 
angle, though the peak is not as pronounced as in cases B2 and B4. This is because the 
OS dileptons can come from either the same gluino (from the leptonic decay of its t and ii 
daughters) or from different gluinos and the cuts tend to favour leptons from top decays to 
those from the decays of stop. 

We show OS lepton event characteristics in Table VII. The SUGRA cases again show 
large jet multiplicity, and large B hadron multiplicity; the latter quantity can again serve 
to distinguish SUGRA events from the MSSM BTW cases. 

Just as with dilepton opening angle, the SUSY source of OS dileptons is correlated with 
the dilepton fiavor asymmetry, defined as 
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^ iV(ee) + iV(/i/x) - Njefi) - iV(e/i) 
^ A^(ee) + iV(/i/i) + Ar(e/i) + Ar(e/i) ^ ' ' 

SS dilepton events are expected to have ~ (replacing ^with i), reflecting uncorrelated 
sources of the leptons; however, OS dileptons may be expected to have a large positive Af 
if many of the OS dileptons come from e.g. Z2 ^ iiZi decay. The asymmetry is listed in 
the last column of Table VII. Indeed, cases Bl and B2 have large Af, while Al and A2 have 
Af consistent with zero. Cases B3 and B4 show intermediate values of Af. Notice that the 
flavour asymmetry qualitatively tracks the azimuthal angle peaking illustrated in Fig. 4. 



D. Same-sign (SS) isolated dilepton events 



The same-sign isolated dilepton signature has been recognized |^JT^ for several years 



now as a particularly clean signal for SUSY, owing to tiny SM backgrounds. It is indicative 
of the Majorana nature of the gluino; in a gluino pair production event, both gluinos may 
decay to charginos of the same sign so that subsequent chargino leptonic decay leads to the 
same-sign dilepton signature. The cross sections after the cuts discussed above are listed 
in Table III, under the SS column. The SUGRA cases yield cross sections between 0.4-3.7 
pb at the LHC, which would yield 4000-37,000 events for an integrated luminosity of lO"^ 
pb~^ per year. Backgrounds are expected to come from ti production, where a lepton from b 
decay is somehow recorded as isolated and from W^W^ production. Our simulation of 300K 
tt(160) events gives a background of 0.01 pb, well below signal whereas the corresponding 
background from same sign WW production has been shown to be ||38|JT6[| ~ 0.035 pb. We 



thus expect that this sample is essentially free of physics backgrounds. 

Each lepton in a SS dilepton event usually originates from the cascade decay of different 
particles produced in the hard scattering subprocess. Thus, the transverse SS dilepton 
opening angle (j){i£') distribution should be maximal for large opening angles. This is shown 
by the distribution plotted in Fig. 4. This behavior, in comparison to that of OS dileptons, 
reflects the source of the dilepton parentage, and hence may be useful in ultimate event 
reconstruction. Furthermore, if SS and OS dileptons had the same source, for instance as in 
gg events, then one would expect roughly equal rates for SS and OS dilepton production. In 
Table III, we present the ratio of OS/SS event rates, where the same cuts have been made 
for each topology, i.e. the cut has been removed for OS dileptons. We see in case 

A2, where dileptons come mainly from t decay in gluino pair events with g ^ tti, that this 
is indeed the case. In case Al, leptonic decays of the top and stop daughters of the same 
gluino are an additional source of OS dileptons. In cases Bl and B2, where OS dileptons 
come frequently come from Z2 iiZi decay, the OS/SS ratio is ~ 6 — 8. Cases B3 and B4 
yield intermediate values of this ratio, in agreement with Fig. 4 and the flavour asymmetry 
values in Table VII. Finally, we show in Fig. 5 the fast and slow lepton px distributions 
for SS events. The hardness of these distributions again reflects the source of the dileptons. 
In particular, the slow lepton distribution in OS dilepton events is typically softer than the 
corresponding distribution in SS events. This is reasonable since there is less energy for 
multi-leptons produced in the same SUSY particle cascade decay. 
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E. Three and Four isolated lepton events 



Complicated superparticle cascade decays can lead to events with 3, 4 and even 5 isolated 
leptons [1^ if gluinos are heavy enough. These cross sections, after the > 100 GeV cut, 
are shown in the 3£ and 4£ columns of Table III. We should warn the reader that in our 
simulation, we have just 10-80 trilepton signal events depending on the case in question and 
just a handful of M events. These numbers could, therefore, contain significant statistical 
fluctuations, and should be viewed in proper perspective. We see that except in case A2, 
~(10'^ — 10^) trilepton events should be expected in a data sample of 10 fh~^ at the LHC. 
However, since gluinos are relatively light within our framework, the M cross sections are 
rather small. SM backgrounds to these multi-lepton topologies are small: for our simulation 
of 300K ti events, no isolated multi-leptons were found, leading to an upper bound on BG 
from this source of < 0.004 ph. Estimates of backgrounds from titi and three and four vector 
boson production also yield rates below the fb level [0]. 

In the M channel, we see a wide range of possible rates for our SUGRA cases. Case A2, 
which is dominated by titi and gg production, where g —>■ tti and ti cZi, produced no 
3£ or 4£ events, corresponding to an upper limit on this rate of ~ 0.02 ph. Case Al, which 
is similar to A2 except now ti bWi, can yield substantial 3i and 4i rates in gluino pair 
events where both tops and stops decay leptonically. Cases Bl and B2, where Z2 iiZi is 
common, can also produce large rates for 3i and 4i events. Cases B3 and B4 do not produce 
large rates for 3£ or 4i events due to a small Z2 £iZi branching ratio (see Table II). 



VI. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have studied the implications of the simplest supergravity grand unified 
model (based on SU(5) gauge symmetry) with radiative electroweak symmetry breaking for 
experiments at the LHC. This necessarily involves assumptions about physics at the ultra- 
high energy scale. Our assumptions, which include the minimality of the particle content at 
and below the GUT scale together with the unification of soft-SUSY breaking parameters 
at the unification scale, are carefully elucidated in Sec. II. In the framework of Minimal 
Supergravity, the SUSY parameters at low energy can all be obtained by renormalization 
group evolution of just four GUT scale parameters together with the top quark Yukawa 
coupling at the unification scale. We have further restricted the range of these parameters 
by incorporating constraints from the non-observation of proton decay [Eq. (2.3)] and from 
the requirement that the relic LSP density is compatible with the present age of our universe 
[Eq. (2.4)]. We have taken at face value the implications of the equality of the 6-quark and 
r-lepton Yukawa couplings at the unification scale (allowing for 10% deviations from exact 
equality due to threshold effects) which, as has been pointed out by several groups, implies 
that the top quark mass is close to its fixed point value. 

As discussed in Sec. HI, these requirements lead to a distinct pattern of sparticle masses 
and mixings. We find, 

• the gluino is light, nig < 400 GeV; 

• the lighter t-squark is substantially lighter than all other squarks, and further, bi is 
lighter than the other squarks; 
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• the Higgsino mass parameter is large, so that the Ughtcr charginos and neutralinos 
are dominantly gauginos, with masses smaller than about 100 GeV; 

• tan (5 is close to unity; 

• the lightest scalar in the Higgs sector is never significantly heavier than 100 GeV. 

We have seen that these mass and mixing patterns lead to characteristic features in 
SUSY events at hadron supercolhders. Of course, these characteristics can be mimicked 
by appropriately adjusting the multitude of parameters of the MSSM. The point, however, 
is that an observation of several of these features could serve as a strong indication of the 
simplest unification scenario. 

Despite the various theoretical and experimental constraints that we have imposed, the 
parameter space is too large to allow a systematic exploration. Instead, we have identified 
six sets of model parameters (that respect all the constraints in Sec. 11) which, we believe, 
capture the diversity of signals that may arise within this framework. The input parameters 
along with the sparticle masses are summarized in Table I. In Cases Al and A2 the gluino 
decays via g tii, while in cases B, gluinos decay via three body decays. ^The main 
difference between Al and A2 is in the stop decay patterns: in case Al, ti bWi while in 
case A2, ti cZy. 

Our main results for the signals from the supergravity framework are presented in Sec. 
V. In order to compare signals from supergravity SU(5) with those of the usually studied 
"minimal models", we have also studied the same signals for the four scenarios labelled 
BTW1-BTW4 defined in Sec. IV, which have a similar value of m-g as in Table 1, but 
have squarks and sleptons exactly degenerate. As usual, we have classified the signals by 
their leptonic content. Table 111 gives an overview of the various signals after all the cuts 
discussed in Sec. V. We see that there are observable signals in essentially all the channels 
shown including the \i and t^l" channels which are usually thought to be swamped by SM 
backgrounds. Unambiguous distinction between all the various supergravity cases, or for 
that matter, even between SUGRA and BTW cases, is not possible, on the basis of cross 
sections alone. 

In order to facihtate this distinction, we have hsted the characteristics of the \i and 
t^t' event samples in Tables V, VI and Vll, respectively. A major difference between our 
supergravity cases Al, A2, B3 and B4, and the BTW cases, is that the supergravity events 
are significantly richer in B hadrons. Of course, whether distinction of scenarios based on 
B identification is actually possible will crucially depend on the ability of the detectors 
to identify displaced vertices. In our calculation, we have assumed a 40% efficiency of 
identifying central B mesons with > 20 GeV in the central region. The reason that 
supergravity events are 5-rich is that one, or more, of the third generation of squarks is 
significantly lighter than the other squarks. In cases Al and A2 where the gluinos decay 
via tti, each gluino pair event contains at least two h quarks from the decays of the top (in 
case Al, the decay of each of the stops also gives 6's). In cases B1-B4, the enrichment of 
the B sample arises in part because 6^ is lighter than the other squarks. As a result, gluino 
and neutralino decays mediated by a virtual 6^ are enhanced. Furthermore, since Z\ and Z2 
are mainly gaugino-like, the Z-mediated decays of Z2 are strongly suppressed. Then, when 
squarks are very heavy, contributions from Higgs mediated Z2 bbZi decays, which further 
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enrich the b sample, become significant. Finally, in cases B2-B4, S's can also come from 
g With decays, although in case B2, the primary h quark is likely to be soft. In cases 
B3 and B4, the primary h quark and, in all three cases, the h quark from top decay should 
have large efficiency to pass our cuts. BTW type events can be very rich in 5's if can 
be produced via cascade decays. As we have mentioned, this is unimportant for the smaller 
gluino masses of interest here (presumably, the kinematic features of the events will give a 
rough measure of gluino mass). We also see from Table V that in cases Al and A2, and also 
in cases B3 and B4, supergravity events can have very high jet multiplicities — about 10% of 
the sample has Uj > 8 whereas this number is somewhat smaller for the BTW cases. For 
cases Al and A2, this is due mainly to the production of top quarks in gluino decays, which 
leads to a large jet multiplicity. Finally, we see that the flavour asymmetry Af defined in Sec. 
V is generally larger in the BTW cases than in the SUGRA cases, except when the gluino 
is very light. This asymmetry, which is a measure of the rate for neutralino production 
in gluino and squark cascades, is greatly reduced when chargino cascades dominate as, for 
instance, in cases Al and A2. 

The features of Tables V-Vll also facilitate some distinction between the various super- 
gravity scenarios. Clearly Al and A2 are readily distinguishable via the high B-multiplicity 
in the Ifj, sample and also from the absence of any flavour asymmetry in the sample. 
The ratio of SS to OS dilepton production rate in Table III may facilitate further distinction. 
Also, the S-multiphcity is significantly higher in Al than in A2. The dilepton cross section 
is considerably higher in case Al (where stops can also yield leptons) than in A2, while 
exactly the opposite is true for the Ifj, signal. Furthermore, since Z2 rarely occurs in the 
Al and A2 cases, case A2 yields no 3£ or 4£ cross section, while Al yields such events only 
from top and stop decay. 

Distinction amongst the various cases B is somewhat more subtle. Bl and B2 are very 
similar, except for the squark mass, as are B3 and B4. A striking difference between (Bl, B2) 
and (B3,B4) is in the flavour asymmetry, in the OS/SS ratio in Table III, in the differential 
distribution shown in Fig. 4, and especially in the large rate for isolated 3£ and 4£ events 
(as has been explained, these are all measures of Z2 production via cascade decays). Some 
distinction may also be possible via the scalar Et distribution since the parent gluino masses 
are signiflcantly different. We should mention that the hard scalar Et cut tends to wash 
out this difference; therefore, this difference may best be studied in the I^j, sample prior to 
this cut. This should be possible since the signal is larger than background even before this 
cut as can be seen from Table IV. The multiplicity of B mesons may provide yet another 
distinction. To distinguish Bl from B2 or B3 from B4 is more difficult. This is not surprising 
since essentially the only difference comes from the squark mass. The clearest distinction 
appears to be in the cross section for the 11 signal and, in the case of B3 and B4, also from 
the If J, cross section. It is encouraging to see that the various scenarios are distinguishable 
even in the messy environment of a hadron collider. 

The reader may have noticed that even though the i-squark is light and even though 
its production constitutes a substantial portion of all SUSY events (see Table II), we have 
not extracted any direct signal for its production. In fact, signals for ti-pair production 
may even be observable at the Tevatron in several of the cases discussed in this paper. 
Since our primary purpose here was to study the distinguishability of various scenarios 
within the SU(5) supergravity framework from one another and from the so-called "minimal 
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models", we used rather hard cuts to obtain clean event samples. The direct signals from 
the lighter sparticles are essentially removed by these cuts and, in our analysis, the existence 
of the lighter stops and sbottom quarks manifested itself in the decay patterns of the parent 
gluinos. It may, of course, well be possible to devise a strategy that is dedicated to look for 
the light scalar top signal at the LHC Charginos and neutralinos may also be searched for 
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In summary, we have examined the phenomenology of the minimal SU(5) supergravity 
grand unified model with radiative breaking of electroweak symmetry. We have delineated 
the range of parameters allowed by constraints from proton decay and from a cosmologically 
acceptable LSP. We have shown that these lead to a characteristic pattern of sparticle masses 
and mixings. We have examined signals from the production of supersymmetric particles 
at the LHC for representative points in the parameter space. We find that these signals 
are distinguishable from the corresponding signals within the usually studied models where 
squarks and sleptons are taken to be degenerate. We further find that by simultaneously 
examining various signals it should be possible to distinguish between the various scenarios 
that are possible even within the restricted framework of supergravity SU(5). Discovery of 
supersymmetric signals at the LHC with patterns corresponding to any of these cases would 
be extremely dramatic, since it would not only signal the discovery of supersymmetry, but 
would provide support for the ideas underlying supergravity grand unification. 

Note added: As we were completing this manuscript, we received a paper addressing some 
of the phenomenological consequences of Yukawa unified "no-scale" models / (^/ . 
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TABLES 



TABLE I. Parameters and masses for six SUGRA cases Al, A2 and B1-B4. 
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TABLE II. (a) Fractions of SUSY particle pairs produced in pp collisions at the LHC; and (b) 
branching fractions of selected decay modes, for six SUGRA cases Al, A2 and B1-B4, where q 
stands for all squarks except stops, and XX stands for all possible chargino and neutralino pairs: 
WiWj, Z,Zj, and WiZj. 
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TABLE III. Cross sections in pb for various event topologies after cuts described in the text, 
for pp collisions at \/i = 14 TcV. The various SUGRA cases arc listed in the first column. The 
OS/SS ratio is computed with the OS dilepton sample before the scalar E-j- cut. 
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TABLE V. Cross sections and event fractions for missing energy plus jets events for various 

SUSY cases at the LHC. 
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TABLE VL Cross sections and event fractions for events containing a single lepton plus jets 
for various SUSY cases at the LHC. 
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TABLE VII. Cross sections and event fractions for events containing opposite sign dileptons 
plus jets for various SUSY cases at the LHC. 
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FIGURES 



FIG. 1. Transverse momentum distribution of fastest jet in multi-jet plus events for four 
SUGRA cases and the sum of backgrounds (histogram), at -y/i = 14 TeV. 

FIG. 2. Distributions for a) scalar Et^ b) lepton plus missing energy transverse mass, and 
c) pT{lepton) in single isolated lepton events for four SUGRA cases Al (dashes), A2 (dot-dash), 
Bl (double-dot-dash) and B4 (triple-dot-dash), along with SM background (solid histogram), at 
■s/s = 14 TeV. Distributions a) and b) are made before scalar Et and mr cuts, while c) is after 
cuts. 

FIG. 3. Distributions for a) scalar Et, b) fast lepton transverse momentum, and c) slow lepton 
transverse momentum in opposite sign isolated dilepton events for four SUGRA cases Al (dashes), 
A2 (dot-dash), B2 (double-dot-dash) and B4 (triple-dot-dash), along with SM background (solid 
histogram), at = 14 TeV. Distribution a) is made before the scalar Et cut, while b) and c) are 
after the scalar Et cut. 

FIG. 4. Distributions for dilepton opening angle in the transverse plane for OS dileptons and 
SS dileptons for SUGRA cases Al, B2 and B4. 

FIG. 5. Distributions for a) fast lepton transverse momentum, and b) slow lepton transverse 
momentum in same sign isolated dilepton events for four SUGRA cases Al (dashes), A2 (dot-dash), 
B2 (double-dot-dash) and B4 (triple-dot-dash), along with SM background (solid histogram), at 

= 14 TeV. 
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